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We show in quite a general way that highly accurate modeling of 
multipath fading responses is possible using low-order complex poly- 
nomials. This applies to all terrestrial radio systems in the chan- 
nelized common carrier bands below 15 GHz, where channel widths 
are 40 MHz or less. The context of the study is a new multipath 
experiment being conducted in New Jersey over a 23-mile path at 11 
GHz. The transmitted signal consists of up to nine tones in a 40-MHz 
bandwidth. These tones are coherently processed, sampled, and dig- 
itized in the receiver and recorded, during fading events, for later 
off-line reductions. Simple routines can be used to determine poly- 
nomial coefficients from these recorded data. This paper describes 
the signal processing and data reduction methods and analyzes them 
to assess the accuracy of polynomial fitting. The analysis uses a 
mean-square error measure and assumes a representative form for 
the underlying response function. Our results predict that the vast 
majority of multipath fading responses can be accurately approxi- 
mated over bandwidths of 40 (62) MHz using first- (second) order 
complex polynomials. 

I. INTRODUCTION 

Multipath fading (hereafter abbreviated mpf) on terrestrial micro- 
wave paths can be a major cause of outage in digital radio systems. 1 " 4 
Numerous efforts have been aimed at understanding, analyzing, and 
correcting this source of disruption, and some have led to new statis- 
tical models for mpf responses. 5-9 

The particular model that inspires the present work approximates 
the mpf response by a low-order complex polynomial in frequency. 8 
For a particular 26-mile path in Georgia, it was shown that a first- 
order polynomial suffices to characterize the fading response in a 25- 
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MHz band centered near 6 GHz. The joint probability distribution for 
the polynomial coefficients was derived for that path, thus permitting 
a complete statistical description of the mpf response. 

Now another experiment is being instrumented, this time for a 23- 
mile path in New Jersey operating in the 11-GHz band. The aim of the 
new experiment is to add to, and in several respects improve upon, the 
data base used to quantify the earlier polynomial model. The improve- 
ments include higher measurement signal to noise ratios (snrs), higher 
sampling rates (20 measurements per second rather than 5), coherent 
processing to obtain phase information (previously absent), and a 
wider measurement bandwidth (40 MHz rather than 25 MHz). 

Given the highly variable nature of multipath fading, such improved 
measurements for a new path in a different frequency band and locale 
should add importantly to our knowledge of this phenomenon. 

The basic design of the experiment can be simply stated: As many 
as nine coherently-phased tones within a 40-MHz bandwidth are 
transmitted from Murray Hill and coherently demodulated in a re- 
ceiver at Crawford Hill; the demodulated tones are sampled, digitized, 
and screened by a desktop computer/controller; and the digitized data, 
if deemed interesting, are recorded on magnetic tape for later off-line 
processing. 

The recorded data will be in a form that facilitates polynomial 
approximation using simple, efficient computer routines. The data will 
be quite general in form, however, i.e., amenable to modeling via any 
mathematical approximation considered promising. 

The present study evaluates the accuracy of polynomial approxi- 
mation, relating it to the experiment parameters and to the methods 
of signal processing and data reduction. Section II describes the signal 
processing in the transmitter and receiver, and derives signal and noise 
relationships used in the subsequent error analyses. Section III de- 
scribes the methods of polynomial fitting to be considered, and defines 
the mean-square error measures that will be used to evaluate them. 

Section IV analyzes the errors in the polynomial fitting caused by 
measurement noise, and Section V analyzes the errors caused by finite 
sampling of the frequency response. In general, the errors increase 
with the bandwidth over which the fitting is done. In analyzing the 
errors caused by finite sampling, we assume a general form for the 
mpf response function that has been applied successfully in other data- 
fitting studies, 5 and assume either worst-case or typical values for the 
function parameters. 

The mean-square error calculations permit predictions of the maxi- 
mum bandwidths for which polynomial fitting is valid. Section VI 
summarizes the results computed, under rather stringent mean-square 
error requirements, for polynomial orders of one, two, four, six, and 
eight. 
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II. SIGNAL PROCESSING ANALYSIS FOR THE MPF EXPERIMENT 

2. 1 Propagation path and radio channels 

Multipath fading responses are to be measured on a 23-mile path 
betwen Murray Hill and Crawford Hill, similar to the one used by 
Crawford and Jakes in their earlier experiments. 10 The transmitting 
antenna at Murray Hill is 655 feet above sea level, and the receiving 
antenna atop Crawford Hill is 425 feet above sea level. An experimental 
license has been obtained to operate over this path in three 40-MHz 
channels within the 11-GHz common carrier band. These channels are 
centered at 11.465, 11.545, and 11.625 GHz. The initial measurements 
will be in the band centered at 11.545 GHz. 

We describe here the signal processing relationships that underlie 
the experiment design. The details of circuitry, components, and 
equipment will be reported separately by those who have developed 
the mpf measurement system. 

2.2 Transmitted signal 

The transmitted signal is created by the two-stage upconversion of 
a baseband signal having the form 

N/2 

b{t) = d + £ d n cos(nAw* + n ), (1) 

n=l 

where N is even and the other parameters will be discussed. The up- 
conversion places the signal in an rf channel centered at radian 
frequency w c = 2w/ c ( f c = 11.545 GHz). Hence, the transmitted signal 
is 

N/2 

V T (t)= £ Jtyn cosiuct + nbut + On), (2) 

n N/2 

where p n is the power of the nth transmitted tone and a total of 
N + 1 tones are transmitted. From (1), we see that p is proportional 
to do and thatp„ (n ^ 0) bears the same proportionality to dl/A. 

The variation of p n with n is clearly symmetrical about n = because 
it derives from amplitude settings of the baseband tones. Nonuniform 
variations of p n can easily be compensated for via baseband adjust- 
ments in the receiver. In Section IV we consider nonuniform variations 
for which receiver noise effects are minimized. 

The frequency spacing between transmitted tones, A/", may be 5, 
10, or 20 MHz. Since the transmission is confined to a channel of 40- 
MHz width, and occupies a bandwidth N&f, we have the constraints 
(N + 1) < 9 tones when A/"= 5 MHz; (N + 1) < 5 tones when A/= 10 
MHz; and (N + 1) < 3 tones when A/= 20 MHz. We will consider the 
four particular combinations N = 2, A/ = 20 MHz; N = 4, A/ = 10 
MHz; N = 6, A/= 5 MHz; and N = 8, Af= 5 MHz. 

Finally, we mention that the N/2 baseband tones are derived from 
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a common 5-MHz reference, and so can be relatively phased in any 
manner desired. For purposes of analysis, we will assume all n 's to be 
zero here; since any phasings in the transmitter are easily undone in 
the receiver, no generality is lost. One criterion for choosing the actual 
dn's is minimization of the peak factor of the rf signal (2). The 
baseband phase adjustments that accomplish this have been derived 
for N = 2, 4, 6, and 8. 11 We will use the resulting minimized peak 
factors in making noise calculations later. 

2.3 Response of the propagation medium 
We denote the complex signal gain of the propagation medium by 
Complex gain at radian frequency co c + <o 



F(«)- 



Gain magnitude (go) during nonfading 



(3) 



The quantity go can be computed from familiar radio path equations. 
Note that co is measured from the center of the channel. During 
nonfading periods, | F(co) | = 1 throughout the channel bandwidth; 
during multipath fading, F(u) varies with co in a randomly time- 
varying manner. 

The function F(co) contains jtwo phase factors of no interest to us. 
One is exp(/0o), where <J>o is the phase shift through the medium at co 
= 0; the other is exp(—jut p ), where t p is the nominal propagation time 
along the path. (For a 25-mile path, t p = 0.13 ms.) The investigation of 
multipath fading can be simplified, with no loss of information, by 
removing these two factors. Thus the function of interest to us is 



//(co) = F(co) exp[y(-«fr> + «*„)]. 



(4) 



The aim of our modeling effort is to find suitable functions for approx- 
imating i/(co), and to statistically characterize the parameters of those 
functions. 

We will see in Section 2.5 that the response function actually 
sampled by the measurement system is 



G(co) =F(co)exp 



->I* + - E 



(5) 



where \p and are the (possibly) random or unknown phases of 
frequency references in the receiver. To obtain samples of the desired 
function [//(co)] from samples of the measured function [G(co)] will 
therefore require performing the operation 



H(u>) = G(co) exp 



j{$ - <;>„) +ju j t p + — 



A* 



AT 



(6) 
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at each of the sampling frequencies [to = 0, ±Ato, • • •, ± (N/2) Ato]. We 
will show later how to accomplish this in the data processing. 

2.4 Decompositions ofF(u), 6(a), and H(u) 

We demonstrate here a useful decomposition for complex response 
functions such as F(w), G(u), and H(u). We will treat only H(w), 
noting that the same mathematics and notation apply to F(u) and 

G(«). 

Since a> is measured from an arbitrary microwave frequency (2-nf c ), 
there is no physical reason to assume complex conjugate symmetry for 
H( to). In its most general form, H(u) can be expressed as 

H(u)=H B (u)+jH b (u), (7) 

where H a (u) and H b (u) are each functions having complex conjugate 
symmetry. Accordingly, we can write 

H a (a) = HaAu) + jH ai (o>), H b (u) = H br (u) + jH bi (o>). (8) 

Real, Imaginary, Real, Imaginary, 

Even Odd Even Odd 

By transmitting and coherently receiving N + 1 tones spaced by Ato, 
one could in theory obtain measurements of the two even functions 
at to = 0, Ato, • • •, (N/2) Ato; and of the two odd functions at to = Ato, 
• • •, (N/2) Ato. [The total number of samples, 2(N + 1), is consistent 
with measuring the amplitudes and phases of the N + 1 received 
tones.] In reality, the receiver obtains these samples for the corre- 
sponding set of G functions, which differ from the H functions if \p and 
9 are not both 0. Obtaining H samples from G samples is discussed in 
Section 2.6. 

Another departure of the receiver outputs from the desired samples 
is the presence of measurement noise. We will defer the introduction 
of noise until Section 2.7. 

2.5 Signal processing in the receiver 
The receiver input at rf is 

N/2 

V R (t)=g t >/2^|F„|cos(to c * + nAto* + <M, (») 

n-N/2 

where | F n \ and <f>„ are the magnitude and phase of F(nAto); go is the 
normal (nonfading) path gain, (3); and we have used (2) with all 8 n 's 
assumed to be zero. 

The signal goes through a two-stage down conversion which amounts 
to quadrature demodulation. That is, two baseband outputs are ob- 
tained which correspond to mixing Vn(t) with 2 coa(w c t + i/O and with 
-2 sin(to c £ + \p). A nonzero value of \p signifies that the rf and if 
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references in the receiver are not in phase synchronism with those in 
the transmitter. 

Each of the baseband signals consists of a dc component plus 
sinusoids at co = Aw, • • • , (N/2)Au. The nth sinusoid in each of these 
signals goes through quadrature demodulation, via the local references 
cos[n(Aut + 8)] and — sin[n(Ato£ + 8)], to produce two more dc 
outputs. These references are all derived from a 5-MHz source in the 
receiver, and nonzero 8 signifies that this source is not in phase 
synchronism with the one in the transmitter. 

Using (5) and ordinary trigonometric identities, the following state- 
ments can be proven: 

(i ) The dc outputs produced by demodulation via 



(2 cos (to c * + *p)\ 
-2 sin(u c t + yp) ) 



are 

7o \ . ( y/2po goG a r(0) 



Qol^UZpogoGbAO)' ""' 



(ii) The dc outputs produced by demodulation via 

2 cos(u c t + «//) 
^—2 sin(w c £ + \f/) 

and cos[/i(Aw£ + 8)] are 

(/-/)„ \ ( y/2p n g Gar(nboi) 
(Q - I)n] - \y/fyngoG br (nbu>) 



n = l, N/2. (11) 

(Hi) The dc outputs produced by demodulation via 



/ 2 cos(u c t + yp)\ 
\-2 sin(w c + \p) ) 



and — sin[n(Aco£ + 8)] are 

'(/- Q) n \ . I V2J^g Gai(nba) 



(Q ~ Q)nJ = \y/2p n goG bi (n*u) ' ' " l ' N/Z ' V? ' ] 

Thus, the dc outputs of the receiver correspond to evenly-spaced 
frequency samples of the four G functions Gor(w), G a i(u), Gbr(u), and 
Gbi(u). Each of the 2(N + 1) dc components is low-pass filtered, with 
a noise bandwidth of 100 Hz; time-sampled 20 times per second and 
quantized with 14-bit precision by a Datel System 256 acquisition 
system;* and passed to a HP 9845A computer for real-time screening. f 



* Datel Corporation. 
f Hewlett-Packard. 
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Any set of samples that seems interesting, or is part of a sequence that 
seems interesting, is recorded on magnetic tape for subsequent off-line 
processing. 

2.6 Relating the G and H functions 

By using (6), the H functions (H ar fa), Haifa), etc.) can be easily 
obtained from the corresponding G functions once A<I> and AT are 
specified. Defining a new function G (to; AO) A Gfa) exp(y'A4>), we first 
perform the matrix operation 



G ar (u>; A0>) 
G 6r (to; AS>) 



cos A<I> 
sin A4> 



-sin AO 
cos A<E> 



G ar (0>) 
G br fa) 



(13) 



An identical equation relates G Q1 (to; A$) and G bi fa; A<D) to G a ifa) and 
G bi fa). The H functions are then obtained from these new G functions, 
for specified A T, as follows: 



Harfa) 
Haifa) 



cos to AT —sin to AT 
sin to AT cos to AT 



G a rfa; A$) 
Gaifa; A$) 



] 



(14) 



An identical equation relates H br fa) and H b ifa) to G br fa; A$) and 

G bi fa; A4>). 

The operation indicated by (13) leads to the result H br (0) = 0, i.e., 
the phase response of the function to be analyzed is forced to zero at 
u = 0. [To see this, combine (6), (12) and (13) for w = 0.] This is a 
welcome simplification in the data and entails no loss of useful infor- 
mation. Fortunately, sin A4> and cos AO are readily obtained from the 
measured G samples, for 

G„r(0) ___ AA _ G ar (0) 



sin A$ = - 



y/GUO) + G'UO) ' 



cos AO = 



VG 2 ar (0) + Gl(0) 



To see this, use (5) at to = and recall that A4» = ^ - </>o.] 

By way of contrast, the value of AT to use in (14) is not so readily 
specified or determined. Yet, to get the full benefit of polynomial 
modeling (i.e., accurate fitting using low-order functions), AT must be 
carefully chosen.* We have arrived at a criterion for choosing AT 
based upon the following data reduction procedure: For a given AT, 
(14) is applied and the resulting H samples are fitted by a finite-order 
complex polynomial inyto. We consider that value of AT to be optimal 
for which the polynomial fitting is best, in some least-squares sense 



* To see why, consider the example G(w) = exp(-ywto), where to a l/B, and B is the 
bandwidth over which #M = G(w) exp(>wAT) is to be fitted. If AT = 0, H (o>) = 
cos ut - j sin uto, and high-order polynomials are needed for accurate fitting over the 
bandwidth. If AT = to, however, //(w) = 1 + >0, and low-order polynomials are quite 
sufficient. 
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defined later. In Section 5.4, we will identify a data-derived measure 
that accurately predicts the optimal A T. 

2.7 Measurement noise 

We have shown that the dc receiver outputs are proportional to 
frequency samples of the function G(co), and that the unwanted phase 
factors that distinguish G(oj) from H(u>) can be removed in the data 
processing. Not so readily removed are the noises associated with the 
digitized outputs. These consist of both additive Gaussian noise from 
the input and components of the receiver and quantizing noise from 
the 14-bit analog-to-digital conversions. 

Receiver noise produces an additive random component for each dc 
output defined by (10) to (12). These 2 (A/ + 1) noises are zero-mean 
and mutually independent. All have the same variance except those 
associated with (10), i.e., n = 0, for which the variance is 3 dB higher. 
These findings follow from the receiver processing described in Section 
2.5. 

We shall now assume tha t eac h of the dc outputs in (10) to (12) is 
adjusted by a factor l/(go V2p„), n = 0, N/2, before being digitized. 
Accordingly, the variance of the Gaussian noise associated with a given 
output sample is 



_ J kTbNp/pogl, n = 0, 



° G \kTbN F /2p n g\ n = l,N/2. (16) 

Table I defines the quantities in (16) and gives values for each. Using 
those data and assuming uniform tone powers, we obtain the following 



Table I — System parameters used in noise analysis 

Parameter Definition and Assumed Value (s) 

g\ Power gain between transmitter output and first receiver amplifier stage; includes 
clear-air path loss and circuit and waveguide losses: 

101og#j = -74dB 
kT Thermal noise density at receiver input: 

10 log kT = -174 dBm/Hz 
b Noise bandwidth of receiver processing for each tone: b = 100 Hz 

Nf Noise figure of first receiver amplifier stage: 

10 log Np = 2 dB 
P p Peak transmitter power, constrained to meet out-of-band emission requirements: 

10 log P p = 14 dBm 
Pf Transmitter peak factor, minimized via phase adjustments of the baseband tones 
before upconversion. For uniform tone powers, 

7.8 dB, N = 2, 
6.3 dB, N = 4, 
6.3 dB, AT = 6, 
6.6 dB, N = 8 
Average transmitter power: P = P p /Pf- For N ■ 8, 
10 log P = 7.4 dBm 
p n Averge power in transmitted tone at frequency ±nAa> from center. 
For uniform tone powers and N = 8, 
101ogp„ = -1.8dBm 



10 log P F = 
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result: ah in dB, for n ^ 0, lies in the range -79 dB ± 2 dB, where the 
precise value depends on N; for n = 0, ah is 3 dB higher. 

Assuming the dc outputs are amplitude-adjusted as indicated, the 
input to every 14-bit A/D conversion is precisely a sample of a G- 
function. During normal propagation, the G samples lie within ± 1 [see 
(3) and (5)]. To provide some room for excess gain, we assume 
quantizer amplitude limits set at ± 1.60 (4-dB margin). As a result, the 
quantizing error for each digitized sample can be characterized as an 
additive noise uniformly distributed on [-A/2, A/2], where A = 2 X 
1.6 X 2~ 14 = 1.95 X 10~\ The quantizing noise variance is then 

a% = A 2 /12 = 3.18 X 10" 9 (-85dB), every sample. (17) 

Comparing this with ah above, we find justification for ignoring quan- 
tization effects. Alternately, they can be accounted for using an ap- 
proximate correction factor given in Section IV. 
To simplify matters further, we introduce the notation 

Har.n A Harinbo)), H br ,n 4 #6r(nA<o), etc. (18) 

These are the quantities produced by (13) and (14) for any specified 
combination of A4> and AT. We account for the noisiness of the H 
samples via the notation 

Har.n 4 Har.n + U.n, H br ,n 4 H br ,n + { br.n, etc., (19) 

where ^ ar ,n, SVn, etc., are the random noise samples. Since the fs are 
produced by phase rotations of the noises associated with the G 
samples, they are all Gaussian, zero-mean, and mutually independent, 
just like the original noises. Moreover, their variances are identical to 
ah, as given by (16). We thus have an accurate, simple description for 
the noisiness of the data to be processed. 

III. POLYNOMIAL FITTING AND ERROR MEASURES 
3. 1 Fitting polynomials to the H samples 

The implicit assumption of the polynomial fitting approach is that, 
over some finite bandwidth B centered on /= 0, the response function 
H(u) can be accurately approximated by a low-order complex poly- 
nomial, i.e., 

M 

H(u>) a H(<o) 4 S (A k +jB k )(ju>) k , | w | < ttB. (20) 

Using (7) and (8), we can break this representation down as follows: 
m m m 

H a (io) = l M/«)*« X MM* + S A k (ju)) k , (21) 

A-0 A-0 A-l 

Even Odd 



Fit to HarM Fit to jH.Au) 
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M M M 

fli(«) = 2 B k {ju) k = 1 MM" + I B k (j u ) k . (22) 

k-0 k-0 k-l x ' 

Even Odd 



Fit to H br {w) Fit toyH w (u) 

The Ak's and Bk's are slowly varying random coefficients; in any given 
measurement interval, they collectively characterize the short-term 
frequency response of the propagation medium. 

The fittings indicated above can be done, for every 50-ms measure- 
ment interval, by using the 2(2V + 1) fl samples obtained in that 
interval. The way the fitting is done depends on the values of M and 
N. We now consider three possible cases. 

Case 1: M = N, with N = 2,4, 6, or 8 

The ^samples obtained using N + 1 tones can be fitted precisely 
using an Nth-order complex polynomial. Thus, when M = N, fitting 
consists of matching each summation in (21) and (22) to the appropri- 
ate H samples at the sample frequencies. The resulting equations for 
the Ak's are as follows (identical equations apply to the Bk's, with 
i/& r 's and Hbi's replacing the Har's and i? a ,'s): 



A k = 



Har.O, k = 0, 

I N/2 

X D%/2,n[H a r,n ~ Har.o], k = 2, 4, • • • , N, 



1 N l 2 
(Aw) , , 



J X 0<A+l>/2.n fiai,n, k = 1, 3, • • • , N - 1, (23) 



where D e i, m and D°, m are the (/, m)th elements of the N/2 X N/2 
matrices [D e ] and [D°], respectively; [D e ] and [D°] are the inverses 
of the matrices [d e ] and [d°], respectively; and the (m, /)th elements 
of [rf e ] and [d°] are 

d' m ,i-(-l) l m a (24) 

and 

d° m ,i=-(-l) l m 2l -\ (25) 

The matrices [D e ] and [D°] for N = 2, 4, 6, and 8 are given in Table 
II. Note for future reference that the derived Ak's and Bk's are weighted 
sums of the H samples. 

Since N can be as high as eight, this method of fitting suggests the 
possibility of eighth-order polynomial modeling. Earlier studies, how- 
ever, suggest that this order is unnecessarily high for bandwidths of 20 
to 40 MHz. 812 Reductions using M = N = 8 may therefore involve 
excessive demands on data storage and analysis, and unduly complicate 
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Table II- 


-Matrices for evaluating eq. (23) 






N 






Matrix 


[D'} 




Matrix [D'] 


2 


[1] 


[1] 


4 






4 
3 

1 

" 3 




l 
Tz 

i 

12 








4 
3 
1 
3 


1 
6 

1 
~ 6 










3 
" 2 




3 
20 




1 
"90 




3 

2 




3 

To 


1 

30 




6 




13 
24 




1 
6 




1 
" 72 




13 
24 




1 

3 


1 
24 








1 
24 




1 
60 




1 
360 




1 
24 




1 
30 


1 

T20 






8 
5 




1 

5 




8 
315 


1 
560 


8 
5 




2 
5 


8 
105 




1 

Tib" 




61 
90 




169 
720 




1 
30 


7 
2880 


61 
90 




169 
360 


1 

To 




7 

720 


8 


29 
360 




13 
360 




1 
"120 


1 

1440 


29 
360 




13 
180 


l 

40 




1 
360 




1 
~ 360 




1 
720 




1 
" 2520 


1 


1 
360 


- 


1 
360 


1 

840 


- 


1 




20160 


5040 



model development and usage. For this reason, we also consider the 
combinations M=N=6,A/=5 MHz; M = N = 4, A/= 10 MHz; and 
M=N=2, A/=20MHz. 



Case 2: M = 2, with N = 8 

Another possibility is to assume a low-order polynomial while using 
all nine tones, in which event the fittings indicated in (21) and (22) are 
done by least-squares methods. Compared with the case M = N = 2, 
this approach requires more data storage and analysis. On the other 
hand, it leads to better fitting accuracy and protects against loss of 
tones caused by equipment problems. Noise effects can be slightly 
worse for this case, despite the noise-averaging it affords, because the 
power per tone must be lower to satisfy the transmitter power con- 
straint. 

The least-squares polynomial fitting method is well known, 1,3 and so 
we give only the results: 



59 


Har.V+2 I fl- 


5* 2 \ 


- 


(26) 


1 4 
Ai -30^L nHa ""' 


(27) 


M " 23U 


iJ 


r 4 

Har.O + 2 £ 
n=l 


(-*> 


ar,n 


(28) 
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The equations for Bo, Bi, and J3 2 are identical, vrithH br , n and# 6 ,,„ 
replacing H ar , n and A a i, n . Again, the A k 's and Bk's are weighted sums 
of the H samples. 

Case 3: M = 1, with N = 8 

The case M = 1 represents the ultimate in modeling simplicity, 
namely, a response function which is first order in Ju. Strong experi- 
mental evidence for such a polynomial has been reported. 8 ' 9 Assuming 
least-squares fitting to nine tones, we find that 



*-{ 



Hor.O + 2 2j H ar ,n 



(29) 



A\ is given by (27); and the same equations give B and Bi, with ■#*,,-.„ 
and ffbi,n replacing ft ar , n and H„i §n . 

3.2 Error measures 

The mpf response function, as reconstructed from the data-derived 
Ak's and Bk's, is 



if 

A-0 



A(u)- 2 (Ak+jBk)(ju) k . (30) 



Since the coefficients are weighted sums of noisy H samples, (19), we 
can write H(u) in the form 

ft(u) = H(u) + Z(u) (31) 

Weighted sum Weighted sum 

of true of random 

H samples noise samples 

Whereas H(u) is noiseless, it is not identical to the true response 
function, H(u), but is a polynomial approximation. 
We now define the error function 

e(w) 4 fi(co) - H(u) 

= [#(«) - #(«)] + Z(w) (32) 

Approximation Noise error, 

error, e(approz) € (noise) 

If B is some bandwidth about f c over which H(u) is to be characterized, 
then the mean-square error for a given polynomial fit can be defined 
as 



1 f B/2 



6 2 4 I J |e( W )| 2 d/ 

i r B/2 i r B/2 

= - \B(u) - H(u) | 2 df+ - \Z(u) | 2 df. 

D J-B/2 ° J-B/2 (33) 

e (approx) c (noise) 
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A useful normalizing quantity for these mean-square errors is the 
mean-square gain of the fading channel, i.e., 

— 1 f B/2 

# 2 A-M \H(<o)\ 2 df. (34) 

B J-B/2 

Radio system experience has shown that, for B > 20 MHz, IP is 
seldom lower than 10 -4 (40-dB total power fade). 

Let us now interpret B as the bandwidth over which the intended 
receiver output (undistorted by mpf) has a roughly uniform power 
spectrum, and outside which the spectral content is small. Accordingly, 
e 2 accurately represents the mean-square error in predicting the output 
signal using #(w) for the mpf response and H 2 accurately represents 
the true mean-square o utput. It is th us r easonable t o say that a given 
polynomial fi t is valid if e 2 (noise) /H 2 ande 2 (approx)/H 5 are both 10" 3 
or less. Since c 2 (noise) is independen t of H(u) , and since IP is seldom 
lower than 10~ 4 , our requirement fore 2 (noise) is that it be 10~ 7 or less. 

IV. MEASUREMENT NOISE ERRORS 
4.1 Method of analysis 

We compute e 2 (noise) by forming Z(u) in (31), using (30) and the 
governing formulas for A k and B k . For simplicity, we use the first-order 
model (Case 3: M = 1, N = 8) to exemplify the approach. 
From (30) and (31), 

i 
Z((o)= S (8A k +j8B h )(ju>) k , (35) 

Ar=0 

where SAo and 5,4 1 are the noises in A Q and Ai, which are computed 
using (27) and (29); and similar definitions and equations apply for 8B 
and 8Bi. There noises are related to the noise components f ar ,o, far.i, 
etc., in the A samples of (27) and (29). Thus, 



M.-J 



4 



far.O + 2 2j !)<>r,n 



SA * = ^T- 2 *Un, (36) 



and similarly for 8B and 8Bi . We now proceed in the obvious way: 
8Ao, SAi, etc., are combined with (35); \ Z(u)\ 2 is formed; and the 
integration in (33) is performed to obtain e 2 (noise). In doing so, we 
make use of the statistical independence among all the noise samples. 
We also use (16) for their variances. The result is of the form 

kTbN F N "S n {B/bf) 



e 2 (noise) = ^V 1 S . (37) 

Ago n-0 Pn 



where, for the case under study, N = 8 and 

4 
81 



4 3 [B\ 2 



n = 0, N/2. (38) 
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For all of the other cases considered, the analytical approach and 
the form of the result are as shown above; only the specific equations 
for S„(B/Af) vary. 

The final step in computing e 2 (noise) is to specify the values of the 
transmitted tone powers. We have considered two approaches, namely 
(£) assuming allp„'s to be equal, adding up to some specified value (P) 
of averag e transm itted power; and (ii) choosing the p„'s so as to 
minimize e 2 (noise), subject to the same average power specification. 
In the first approach, we apportion power according to the rule 

P»-Jf+l> aUn. (39) 

In the second approach, we use the method of Lagrangian multipliers 
to minimize (37), subject to the constraint 

N/2 

po + 2 I p n = P. (40) 

The result in that case is 



J(l + 8 n )S n (B/Af)P 

Pn = , = , all n, (41) 

V2S (B/A/) + 2 £»3 y/Sn(B/lf) 

where 8„ is the Kronecker delta function. Note that the optimal 
variation of p„ with n depends on the bandwidth ratio, B/ A/. 

We have reduced (37) to numerical results using the parameter 
values in Table I. We now present our findings. 

4.2 Results for Case 1{M=N=2,4,6,or8) 



Figure 1 shows curves of e 2 (noise) vs. B for the various possibilities 
under Case 1. We can make the following observations: 

(i) The noise errors and their rate-of-growth with B depend strongly 
on the polynomial order M , both being less for lower M. From a noise 
standpoint then, M should be chosen to be low. The competing factor 
influencing this choice is the approximation error, which we consider 
later. 

(ii) Optimizing p n does not improve significantly on the use of 
uniform tone powers. There is no compelling reason, therefore, to 
taper the p n 's. On the other hand, doing so entails no price in com- 
plexity and may offer benefits, e.g., tapering of p n may reduce out-of- 
band spurious tones caused by transmitter nonlinearities. We will not 
explore this topic here, but present in Table III some examples of 
optimal p n variat ions. 

(Hi) Values of e 2 (noise) less than 10~ 7 can be attained, using either 
uniform or optimal p„'s, for bandwidths up to 34 MHz or more, 
depending on M. For purposes of modeling, errors of this magnitude 
can be regarded as negligible, as noted in Section 3.2. 
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UNIFORM TONE POWERS 

OPTIMAL TONE POWERS 



5 -70 




30 



40 50 60 

BANDWIDTH, B IN MEGAHERTZ 



70 



Fig. 1— Mean-square noise errors, e 1 (noise) in dB, for Case 1 (M = N = 2, 4, 6 or 8). 
For each N, the curve for optimal tone powers (dotted) merges with the one for uniform 
tone powers (solid). 



4.3 Results for Cases 2 and 3{M = 2, N = 8 and M = 1, N = 8) 



Figure 2 shows curves of e 2 (noise) vs. B for Cases 2 and 3. These 
results show the effects of noise when first- and second-order polyno- 
mials are formed via least-squares fitting to nine received tones. For 
purposes of comparison, Fig. 2 also repeats the results for M = N = 2 
under Case 1, i.e., second-order matching to three received tones. We 
note the following: 

(i) As before, fitting H(u) with a lower-order polynomial leads to 
smaller noise-related errors. 

(ii) Comparing the two ways of getting M = 2, noise errors are 



Table III— Optimal tone powers for Case 1 , B = 40 MHz (in dB 
above power of central tone) 



Frequency 



±5 MHz 



±10 MHz 



±15 MHz 



±20 MHz 



-0.80 
-0.36 



-0.10 
-3.30 

-1.36 



-7.60 
-2.63 



-3.02 
-4.05 



-7.95 
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-80 



-85 - 



20 



30 



40 50 60 

BANDWIDTH, B IN MEGAHERTZ 



Fig. 2 — Mean-square noise errors, e 2 (noise) in dB, for Cases 2 and 3 (M = 2, N = 8 
and M = 1, N = 8). Results shown are for uniform tone powers; results for optimal tone 
powers are lower by less than 0.5 dB. 



sometimes smaller when N = 2 and sometimes smaller when N = 8. 
The explanation lies in two offsetting factors: Using nine tones leads 
to less power per tone, but affords averaging benefits not available 
when using three tones. Depending on B, one or the other of these 
factors dominates. 

(Hi) F or all three low-order polynomial approaches considered, 
e 2 (noise) is less than 10" 7 for bandwidths up to 90 MHz or more. 

4.4 Generalizations and extensions 



We emphasize that e 2 (noise) in (37) scales readily in the quantities 
bNp/Pgo and B/Af. Therefore, the curves of Figs. 1 and 2, however 
particularized to the present experiment parameters, can easily be 
scaled to reflect other conditions. 

With this in mind, we can account in a simple and fairly accurate 
way for both quantization noise (previously neglected) and changes in 
P, go, b, and iV>. We assume the quantization and Gaussian noise 
powers to be additive, invoke (16) and (17), and obtain the following 
rule: If bN F /P gl is chan ged by a factor C from the one obtained using 
Table I, then e 2 (noise) should be adjusted by the factor (C + 0.25). 
For C = 1 (no change in bNp/Pgo), quantizing noise adds roughly 25 
percent (1 dB) to the mean-square noise error. 
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V. APPROXIMATION ERRORS 
5.1 Method of analysis 

The approximation error defined in (32) results from sampling H(u) 
at a finite number of discretely-spaced frequencies. The mean-square 
error, as defined in (33), depends on B, A/", M, N, and— unlike the noise 
error— on the specific function being approximated. We consider a 
highly accurate polynomial approximation to be one for which 77, 
defined by 

1, A e 2 (approx)/*P, (42) 

is less than or equal to 10" 3 . Our aim here is to calculate tj using a 
function for H (w) that is "representative," in a sense to be described 
shortly. The method of computation, given H(u), follows directly from 
the relationships in Section III. 

The basic form assumed for H(u) is the one used by W. D. Rummler 
in previous studies. 5,69 To fit measured amplitude responses in a 25.3- 
MHz bandwidth near 6 GHz, Rummler expressed H(u) by the function 
a[l - b exp{j0) exp(-y'wAT)], where a, b, B, and At are function 
parameters to be chosen. For eaclLof 24,920 response measurements, 
these four parameters were chosen to give a least-squares functional 
fit to the data. Error analyses of the results (Fig. 17 of Ref. 5) revealed 
a high degree of fitting accuracy over nearly all of the data records. 
For this reason, we use the same function here to represent the actual 
underlying response of multipath fading. 

Rummler also analyzed the effect of fixing the delay parameter At 
at 6.31 ns, and choosing (a, b, 6) to fit each data record. While not as 
good in all cases, the "fixed delay" model was found to be highly 
accurate over at least 98 percent of the data base. With this in mind, 
we conjecture that the above function, with At < 10 ns, can be used to 
represent the vast majority of mpf responses. We will invoke this 
conjecture later in tabulating bandwidths over which polynomials of 
specified order can be used to model multipath fading. 

Our approach in what follows will be to find tj for specified B and 
At, with b and 6 chosen to give worst-case polynomial fitting results. 
This should provide conservative estimates of the limitations of poly- 
nomial fitting. We will also assume a phase adjustment to the function 
used by Rummler, to enforce consistency with the data processing 
approach described earlier. In particular, we specify the response 
function to be 

H(u) - o(l - be je e- j ^ T ) exp[y(AO + «AT)], (43) 

where A<I> and A T were defined in Section II. As noted there, the data 
reduction process chooses AO in such a way that H (0) is made real, 
and chooses AT in such a way that tj is minimized for a given Mth- 
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order p olynomial fit. Acc ordingly, we specify that exp(y'AO) = (1 — 
be~ je )/\ll + b 2 — 2b cos0; search over AT so as to minimize -q for a 
given fitting polynomial and given (b, 0); and search over b and so as 
to maximize that minimized t/.* Thus, we obtain the measure 



rjo k Max { Min(Tj) > , worst-case tj, given At and B. (44) 
(M) [ ir J 

We have done this for each of the different combinations of M and N 
described earlier, using B At as a variable and B as a parameter. We 
now present our findings. 

5.2 Results for Case 1{M=N=2,4,6,or8) 

The worst-case combination of (b, 0) for any B, At, and polynomial 
order is found to be (1.0, 0.0). This combination corresponds to total 
fading at the band center (to = 0). Also, the optimal AT for b = 1.0 is 
precisely At/2. 

The variations of tjo with BAt, with M and B as parameters, are 
shown in Fig. 3. As expected, t/o is a strongly decreasing function of M , 
a strongly increasing function of BAt, and a weakly increasing function 
of B alone for given BAt. 

A useful empirical formula for these results, accurate to within 40 
percent for B < 80 MHz and BAt < 1.2, is 

t,o = k(M, B)10- m (BAt) 2M [1 + «(BAt)], (45) 

where a = 0.57 and k(M, B) is a fairly mild function of M and B. 
Numerical values for k (M, B ) are given in Table IV. 

5.3 Cases 2 and 3 (M = 2, N = 8 and M= 1, N = 8) 

Curves of tjo for these cases are shown in Fig. 4, where the results for 
M = N = 2 are repeated for comparison purposes. There i6 an evident 
improvement when second-order polynomials are derived using nine 
tones instead of three. An empirical formula for these results (again, 
accurate to within 40 percent for B < 80 MHz and BAt < 1.2) is given 
by (45), where a = 0.32; the power of (BAt) is replaced by 4; and 
numerical values for k (M, B ) are given in Table IV. 

5.4 Optimizing A T using measured data 

In minimizing t/ with respect to AT, (44), we used computer search 
procedures and exploited our assumed knowledge of the underlying 
H( to). In an actual measurement the latter is not possible, as the only 
information available for optimizing AT is the set of measured G 
samples. 



* The amplitude factor a has no impact on fitting accuracy and so is set to unity for 
purposes of this study. 
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Fig. 3— Relative mean-square approximation errors, tj in dB, for Case \(M= N =2, 
4, 6, or 8). The variable is BAt, where At is the delay spread quantity in Eq. (43). 
Bandwidth, B, and M are parameters. 



By experimenting with different strategies, we have identified a 
computationally efficient procedure that uses these samples to opti- 
mize AT. It consists of computing C 2 A (Al + Bl) as a function of AT 
and choosing that A T for which C 2 is a minimum. The procedure is 
even simpler than it might seem; for, if A n , is A n computed for AT = 
0, then A 2 for any other A T is just 

M = A 2 ,o + Ai,oAT + fcAcoAT 2 , (46) 

and an identical relationship applies to B 2 . These results can be 
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Table IV— Values of k(M, B) in eq. (45) 



Case 


Case 1 
(M = N) 


Cases 2 and 3 
(N = 8) 


b\m 


8 


6 


4 


2 


2 1 


40 MHz 


1.31 x 10" 4 


0.138 


3.32 


6.39 


1.47 1.51 


60 MHz 


1.95 x 10~ 2 


0.734 


5.31 


6.98 


2.56 1.78 


80 MHz 


5.83 X 10~ 2 


1.14 


6.17 


7.20 


4.26 2.26 



derived using (6), (21), and (22), and the power series expansion of 
expOcoAT). 

Thus, "optimizing" A T amounts to rinding the minimum of a fourth- 
order polynomial in that quantity. This can be done very efficiently 



O -20 - 



5 -30 - 



-60 




Fig. 4 — Relative mean-square approximation error, tjo in dB, for Cases 2 and 3 
(M = 2, N = 8 and M = 1, N = 8). Same variable and parameters as in Fig. 3. 
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using computer algorithms. Our investigation of this procedure con- 
sisted of minimizing C 2 with respect to A T for numerous (6, 0, At), and 
using the resulting values of A T to compute tj. We obtained remarkably 
similar answers, in most cases, to those obtained with A T rigorously 
optimized [i.e., via computer search using the known function tf(<*>)]. 
The worst departures occur when tj is either so large as to be of no 
interest or so small that the error increases do not matter. 

Other measures derivable from the data, such as [(A 2 —A?) 2 + (B 2 
— Bi) 2 ], may be even more appropriate Quantities such as this, or Co 
as defined above, are reliable indicators of the "curvature" in /f(w); 
choosing M to minimize them should therefore maximize the fitting 
accuracy of low-order polynomials. Our results give empirical support 
for this principle. 

VI. SUMMARY AND CONCLUSION 

The major findings of this study are summarized in Table V. The 
first row gives, for each combination of M and N conside red, the 
bandwidth below which polynomial fitting yields e 2 (noise) < 1(T 7 . 
These bandwidths are derived for uniform tone powers and the param- 
eter values listed in Table I. They cannot be increased much without 
dramatic (and unlikely) improvements in transmitter power, noise 
bandwidths, and quantizing precision. 

The second row of Table V gives the corresponding bandwidths 
below which tj < 10~ 3 . They are derived for the response function (43), 
with b and having worst-case values and At = 10 ns. The third row 
gives the maximum bandwidth which satisfy both mean-square error 
requirements. Since these requirements are quite stringent, we present 
in the fourth row the results of relaxing both of them by 6 dB (factor 
of four). 

These tabulations suggest that the vast majority of mpf responses 
can be accurately approximated, over bandwidths of 40 (62) MHz, by 

Table v — Bandwidths, in MHz, over which H(w) can be accurately 

characterized 



Case 


Case 1 
(M = N) 


Cases 2 
(N- 


and 3 

= 8) 


Requirement^, M 


8 


6 




4 


2 


2 


1 


e 2 (noise) < 10"' 


40 


34 




55 


96 


90 


180 


tj < 10- ' 


140 


125 




96 


35 


48 


28 


Both requirements 


40 


34 




55 


35 


48 


28 


Requirements 
relaxed 6 dB 


43 


37 




62 


48 


62 


40 



Notes: 

(1) Results assume At = 10 ns. 

(2) Results in excess of 80 MHz are extrapolated estimates. 
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complex polynomial functions of just first (second) order. We surmise, 
then, that the appropriate modeling approach is to form first- and 
second-order polynomials using nine-tone measurements {N = 8) and 
least-squares data fitting. With these low orders, there is maximal 
simplicity in developing and using the model, with no substantial 
sacrifice in accuracy. By using nine tones, moreover, polynomials of 
higher order (up to eight) can also be deduced from the data, thereby 
enabling the theoretical predictions given here to be fully tested. 
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